Nanosecond quantum state detection in a current biased dc SQUID 
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This article presents our procedure to measure the quantum state of a dc SQUID within a few 
nanoseconds, using an adiabatic dc flux pulse. Detection of the ground state is governed by standard 
macroscopic quantum theory (MQT), with a small correction due to residual noise in the bias 
current. In the two level limit, where the SQUID constitutes a phase qubit, an observed contrast 
of 0.54 indicates a significant loss in contrast compared to the MQT prediction. It is attributed to 
spurious depolarization (loss of excited state occupancy) during the leading edge of the adiabatic 
flux measurement pulse. We give a simple phenomenological relaxation model which is able to 
predict the observed contrast of multilevel Rabi oscillations for various microwave amplitudes. 




PACS numbers; Valid PACS appear here 



I. INTRODUCTION 

The last few years have shown significant breakthrough 
in the domain of superconducting Josephson qubits. Co- 
herent manipulations have been achieved on a single 
qubit with various designs [l|, [3, H, 0, Q • Recently, cou- 
pling between two qubits has been demonstrated 
and quantum logical gates have been realized [1, [iC 

The dc SQUID is a versatile quantum circuit display- 
ing several energy levels with tuneable anharmonicity. 
In the high anharmonicity limit, it becomes a phase 
qubit whose decoherence has been characterized 
With smaller anharmonicity, microwave excitation in- 
duces complex dynamics which involve many levels [T^ . 

Readout of the quantum state is one of the key issues 
in quantum circuits operation. As proposed in [l^ and 
experimentally implemented in our measurement of 
the quantum state of the dc SQUID is performed within 
a few nanoseconds, using adiabatic dc flux pulses. De- 
pending on the SQUID state, the pulse implements a 
selective switch to the resistive state. The dc SQUID 
design thus combines two functionalities: the quantum 
circuit and the state measurement. A fast measurement 
allows for a time domain study of fast dynamics. A sim- 
ilar measurement procedure for a phase qubit, developed 
independently, has been used to investigate the coupling 
to spurious two level fluctuators [14 1 and more recently 
the coupling between two qubits [8|. Moreover, speed 
also minimizes contrast losses due to depolarization dur- 
ing the measurement stage. In this article, we present 
a complete characterization of our measurement perfor- 
mance, both in the two level and multilevel limit. 

The paper is organized as follows. Section I intro- 
duces the physics of the SQUID and discusses its op- 
eration as a quantum circuit. The experimental imple- 
mentation is presented in Section II, focusing on the mi- 
crowave part which is a key point of the experimental 



setup. In the third part, results on the ground state de- 
tection are presented. Excepting a small correction due 
to residual noise on bias parameters, the experiment is 
consistent with macroscopic quantum tunnelling theory. 
Detection in the two level limit and measurement in the 
multilevel limit are then successively presented. A simple 
relaxation model successfully accounts for measurement 
performances. To finish, section IV contains concluding 
remarks and improvement perspectives. 



II. THE CURRENT BIAS DC SQUID 

The dc SQUID, presented in Fig. [Ija), consists of two 
identical Josephson junctions (JJ), each with a critical 
current Iq and a capacitance Cq. The junctions are em- 
bedded in a superconducting loop of total inductance Lg 
which is unequally distributed between the two arms of 
the SQUID. The asymmetry in inductance is character- 
ized by the parameter rj — (Li — L2)/Ls, where Li and 
ii are respectively the inductance of the first and sec- 
ond arm. In our circuit, the Josephson energy Ej = 
{^o/2tt)Iq is more than 3 orders of magnitude greater 
than the Cooper pair Coulomb energy Ec ~ (2e)^/2Co 
(<f>o — h/2e is the superconducting flux quantum). In 
this limit, the two superconducting phase differences ipi 
and ip2 across the two Josephson junctions are the nat- 
ural variables to describe the dynamics of the system. 
This two-dimensional phase dynamics can be treated as 
that of a fictitious particle of mass m = 2Co($o/27r)^ 
moving in a two-dimensional potential [isj 



U{x,y) — 2£'j[— (j^)a;— cosxcosy— I 
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shape of the potential, which contains valleys and moun- 
tains, is experimentally controlled through the bias cur- 
rent /fc and the bias flux threading the supercon- 
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FIG. 1: (a) Electrical schematic of the dc SQUID. In the 
superconducting state, this circuit can be seen as two induc- 
tively coupled resonators. The two pure Josephson dipoles, 
symbolised by crosses, introduce two non-linear inductances 
controlled by the bias current h and the loop current which 
tends to screen the external applied flux (b) Contour 

plot of the potential surface U{x,y) versus x — ^{'■Pi + 1^2) 
and y = ^{fi — f2), for a typical bias point; (x) and [y) 
scales are identical, energy scale is arbitrary. Along the es- 
cape direction [tp) , which joins the local minimum M and the 
only saddle point S, the potential is cubic (c). The quantum 
states under study are vibration states within the potential 
well. In the transverse direction (ipx), the potential is roughly 
quadratic with a bottom frequency which is about twice Vp. 
Coupling between the escape and the transverse direction is 
neglected. 



ducting loop. Here, the inductance coupling parameter 
b — $o/(27rLs/o) is of order unity. With the values of 
/f, and <I>t, used in our experiments, the potential sur- 
face displays only one type of local minimum. They are 
separated from each other by a saddle point, where the 
particle can escape A contour plot of the poten- 

tial surface for a typical bias point is presented in Fig. 
[Ijb). Along the escape direction, which makes an an- 
gle 9 with the x axis, the potential is cubic and is com- 
pletely characterized by its bottom frequency Vp and a 
barrier height AC/ [Fig. HJc)]. These two quantities de- 
pend on the magnetic flux and vanish at the SQUID 's 
critical current /c We assume complete separation 
of the variables along the escape direction and the trans- 
verse one by neglecting the coupling terms between these 
two directions. In this approximation, the phase dynam- 
ics of the SQUID along the escape direction is similar to 
the phase dynamics of a current-biased single Josephson 
junction. The parameters AC/ and I'p, explicitly given 
in the appendix 1, are renormalized, thereby taking into 
account the two-dimensional nature of the potential. 

For bias currents h < Ic, the particle is trapped in a 
local minimum and its quantum dynamics is described 
by the Hamiltonian 

Ho = ^hi^p [pl + if^] - ahvpif^. (2) 



(p = {2-K mvp / h)4i represents the reduced phase oper- 
ator, associated to the phase ^ = cos^i -|- suvOij along 
the escape direction, is the reduced momentum con- 
jugated to If. The relative magnitude tr of the cubic term 
compared to the harmonic term is related to the barrier 
height AC/ = hvp/^Aa^. For /{, below Ic, several low- 
lying quantum states are found near the local minimum. 
These states, describing the oscillatory motion within the 
anharmonic potential, are denoted |n) for the nth level, 
with n = 0, 1,2, ... The corresponding energies E^i were 
calculated in Ref. The Bohr frequency associated 

with the transition n ^ k\s denoted Vnk in the following. 
Since energy levels are trapped in a finite height potential 
barrier, they are metastable. Escape out of the potential 
well is possible both by tunnelling and thermal activa- 
tion. For temperatures well below T* = hvp/ {2TrkB), the 
so-called cross-over temperature, the thermal process can 
be neglected [l^. On the timescale of the experiment, 
only the three highest energy levels, closest to the top of 
the barrier, undergo a significant tunnelling effect. Lower 
energy levels are stable. 

Our procedure to perform quantum experiments con- 
sists of the repetition of an elementary sequence, which 
is decomposed into four successive steps. A bias cur- 
rent If, is first switched on through the SQUID at fixed 
magnetic fiux $6. The working point {h, $b) defines the 
geometry of the potential well. The Jf, risetime is slow 
enough to induce adiabatic transformation of the poten- 
tial and the circuit is initially in the ground state |0). 
A microwave flux pulse characterized by its frequency v 
and amplitude manipulates the quantum state of the 
system. It induces a time-dependent perturbation term 
— /la^uiV^ cos(27ri^i)(^ in the Hamiltonian which couples 
neighbouring levels. In this expression, a^u; is directly 
proportional to and corresponds to the frequency of 
the Rabi oscillation between |0) and |1) for an excitation 
tuned on vqi (see appendix 2). Then, a dc flux pulse of 
amplitude is applied and brings the system to the 
measuring point {Ib,^b + ^m)- This flux pulse adiabat- 
ically reduces the barrier height and allows tunnelling 
escape of the localised states to finite voltage states dur- 
ing a time At of the order of a few nanoseconds. As 
discussed in Ref. with precise adjustment of $,„ 

it is theoretically possible to induce a selective escape of 
excited states. Because the SQUID is hysteretic, the zero 
and finite voltage states are stable and the result of the 
measurement can be readout by monitoring the voltage 
Vs across the dc SQUID. Escape out of the well corre- 
sponds to the detection of a voltage which is twice the 
superconducting gap of the circuit material. The current 
bias lb must be switched off to reset the circuit in the 
ground state. This elementary sequence is repeated to 
extract the occupancies of the excited states. 

This article is devoted to the characterisation of the 
nanosecond measurement procedure. The experiments 
presented here only involve levels |0) to |3). Denoting 
Pn the population of the nth level, we thus restrict dis- 
cussion to states characterized by the set of occupancies 
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{po,Pi,P2,P3}', generalization to higher energy levels is 
straightforward. For such a state, the escape probability 
Pe out of the potential well reads: 

Pe=P^ + j2iP"-Pe)^Pn- (3) 

Pei^m) = 1 — exp[— r„(<i>m)Ai] is the escape probabihty 
of the state \n). From the measurement point of view, 
P" correponds to the probability to detect this level when 
it is fully populated. The tunnelling escape rate F^ de- 
pends on the geometry of the potential barrier which is in 
turn fixed by the measuring point (/ft, + Escape 
rates are obtained from semiclassical calculations. Par- 
ticularly, Fq is given by the usual expression for macro- 
scopic quantum tunnelling of the ground state of a single 
Josephson junction, 

ro = 12V^-pexp(-fA^), (4) 

and the contrast in escape rate between two successive 
levels is F„+i/F„ — exp(27r) « 535 ^1§\. In the following, 
we refer to the plot of P" versus as the escape curve 
of level |n). 

In the two level limit, the detection contrast between 
states |0) and |1) is P^ — P" . The measuring amplitude 
which optimizes the detection probability allows a detec- 
tion contrast of 0.98 and corresponds to P" = 0.02. With 
these settings, in the general case where more than two 
levels are populated in the system, the escape probability 
is practically equal to the population of excited states: 

Pe == 0.02 + 0.98pi+p2+P3- (5) 

Moreover, it is possible to find a lower pulse amplitude for 
which the ground and first excited states remain trapped, 
but states |2) and |3) escape with close to 100% probabil- 
ity. Repeating this procedure with successively smaller 
pulse amplitudes, one can measure pi + P2 +P3, P2 +P3 
and p3 and therefore reconstruct all levels occupancies. 

Our measurement is one of the fastest implemented in 
Josephson quantum circuits. The rise time of the mea- 
suring pulse is 1.6 ns and the duration Ai can be as short 
as 1.5 ns. This speed allows a time domain study of fast 
relaxing phenomena and tends to limit the loss of mea- 
surement contrast due to relaxation. The question of the 
possibility of nonadiabatic effects, namely transitions be- 
tween energy levels which modify the levels occupancies, 
has been recently addressed for a rf SQUID which is a 
very similar device [13] ■ Calculations show that the adi- 
abatic error is negligible for the pulses used in our exper- 
iments: it is less than 10"'^ for At = 1.5 ns and decreases 
very rapidly for slower pulses. 

In terms of speed, escape to the running state of a tilted 
washboard potential may also have an advantage over 
the method used for the rf SQUID 1^. In that case, the 
fictitious particle escapes to another trapped state, which 
can interact strongly with the qubit levels if relaxation 




FIG. 2: Principle of the SQUID operation, (a) Schematic 
of the electrical circuit. The capacitively shunted JJs are 
symbolised by a square with diagonals. In our experiments, lb 
is first switched on through the SQUID at fixed magnetic flux 
"!>[,, defining the shape of the potential well. A high frequency 
flux signal is then applied on the superconducting loop with a 
microwave antenna. A digital sampling oscilloscope record of 
this signal is presented in (b). It is composed of a microwave 
pulse which induces transition between adjacent energy levels, 
followed by a measuring dc pulse which adiabatically reduces 
the barrier height and induces a selective tunnelling escape 
of excited states. Escape is detected by a voltage Vs across 
the circuit, amplified at room-temperature. At the end of 
the sequence, h is switched off and the circuit is reset to the 
superconducting ground state. 

is insufRcient. This leads to retraping of the particle in 
the qubit states. This has been investigated numerically 
in Ref. [l^]. Assuming an infinite Ti, the retrapping 
probability was found to be greater that 15 %. In our 
case, the qubit excited state escapes to an unbound state 
of the tilted washboard. Retrapping is therefore much 
less likely, and faster measurement times possible. 

III. EXPERIMENTAL SETUP 

The measured circuit was realized using e-beam lithog- 
raphy and shadow evaporation of aluminium on an oxi- 
dized silicon chip. A SEM picture of the circuit is pre- 
sented in Fig. ^ih). The SQUID consists of two large JJs 
of 15 /im^ area enclosing a 350/xm^-area superconducting 
loop. Each Al/AlOx/Al junction has a critical current 
/o = 1.242 /iA and a capacitance Co = 0.56 pF. The 
two SQUID branches combine to give a total loop induc- 
tance Ls — 280pII with asymmetry parameter rj = 0.414. 
These electrical parameters are obtained by combining 
the analysis of two independent measurements. First, 
Ls, rj and an approximate value of Cq are extracted 
from macroscopic quantum tunnelling measurements, as 
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discussed in Ref. [H, [T^. Second, spectroscopic mea- 
surements of the |0) |1) transition frequency uqi as 
a function of h are used to determine a precise value 
for Co [m [13] ■ AH measurements are performed with 
the sample at a temperature of 30 mK in a dilution 
fridge. Thermal energy is then small compared to the 
relevant energy scales of the circuit, namely the Joseph- 
son energy Ej = ks x 29.6 K and the vibration quantum 
\/EjEc = fcs X 442mK, with ks the Boltzmann constant. 

To avoid spurious microwave resonances, the chip is 
mounted in a shielded copper cavity whose cut-off fre- 
quency is above 20 GHz. The immediate electromagnetic 
environment of the SQUID, located in the microwave cav- 
ity, is designed to decouple the circuit from the external 
room temperature classical electrical apparatus. It con- 
sists of two cascaded filters. The first and most impor- 
tant includes a 9 nH on-chip inductance introduced by 
long on-chip superconducting thin wires, and a gold thin 
film 150 pF parallel capacitor. The second filter consists 
of the bonding wires, with an estimated inductance 3nH, 
and a surface mounted 2nF capacitor. This electrical en- 
vironment displays two low frequency resonances whose 
measured frequencies are consistent with the above pa- 
rameters [2l|. As illustrated on Fig. E^a), special care 
was taken in bias line filtering, and the experimental 
setup includes several low temperature filters. Bias flux 
is applied with a copper coil cooled down to 30 mK. Ex- 
ternal low frequency flux fluctuations are screened with 
soft iron and mu-metal at room temperature and a super- 
conducting lead shield at 1.5 K. The room temperature 
high frequency signal is guided by 50 Q. coax lines and 
attenuated twice by 20 dB at 1.5 K and 30 mK, before 
reaching the SQUID through an on-chip short circuit. 
The antenna is located 10//m from the SQUID, leading 
to a 1.5 pH mutual inductance coupling. The /iw line is 
terminated by an inductance estimated to be 2 nH which 
introduces a supplementary 6 dB attenuation at 10 Ghz. 
Special care was taken to minimize the coupling between 
the high frequency antenna and the bonding wire loops 
in the direct electrical environment of the SQUID. 

The high frequency signal used in our experiments re- 
sults from the combination of a microwave excitation 
pulse and a measurement dc pulse. The /zw pulse is gen- 
erated by mixing continuous microwaves with a shaping 
dc pulse. We define the risetime as the duration to reach 
95 % of the peak amplitude starting from 5 %. The dura- 
tion is defined as the time interval where the pulse am- 
plitude is at least 95 % of the maximal amplitude. With 
this convention, the shaping pulse has a 0.8 ns risetime 
and a duration tuneable from 1 ns to several millisec- 
onds. The measurement pulse has a 1.6 ns risetime and 
a minimum duration of 1.6 ns. Both microwave shaping 
and measurement pulse are produced by two outputs of 
one dc pulse generator, allowing precise control of tim- 
ing. In Fig. [U^c), a digital sampling oscilloscope record 
of the signal is presented. As a demonstration of signal 
generation performance, the duration of the microwave 
and measurement pulse have been reduced to the mini- 
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FIG. 3: (a) Schematics of the experimental setup with tem- 
peratures indicated on the right, 
order low-pass rlc filter. 



F-rlc 



stands for a second 



F-th 



is a 1 m lossy thermocoax and 
F-vr I is composed of two cascaded tt filters and two meters 

is a 50 f2 attenuator, (b) Scan- 



-20 dB 



of lossy thermocoax. 
ning electron microscope picture of the SQUID. The circuit 
is obtained by shadow evaporation of two 15 nm aluminium 
layers; Josephson junctions appear in light grey. The extrem- 
ity of the 500 /im long thin current bias line is visible on the 
left: their width is between 200 nm and 300 nm. The end of 
the on-chip pw antenna appears on the right. 



mum; the whole sequence is then finished in less than 
7 ns. Excepting specific measurements, the delay be- 
tween microwave and measurement is kept as short as 

1 ns to minimize depolarization due to population relax- 
ation. The microwave amplitude a^^, can be tuned from 
to 300 MHz. Depending on the bias point, the measur- 
ing pulse amplitude ranges between a few 10~^ $0 find 

2 X 10^1 $0- 

To limit heating in the SQUID voltage state, a specific 
electronic circuit cuts the bias current to zero as soon 
as a finite voltage across the SQUID is detected. The 
voltage readout time, about 10 /iS, is determined by the 
bandwidth of the strongly filtered voltage lines. With 
a wait of 1 ms after each voltage measurement to allow 
the sample to cool back down to base temperature, the 
overall repetition rate is about 1 kHz. The experiment 
sequence is repeated up to 5000 times such that the sta- 
tistical noise on has a theoretical standard deviation 
lower than 4 x 10"'^. 
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FIG. 4: Detection probability of the ground state versus the 
amplitude of a 4 ns measuring pulse. The fit of experimental 
data (o) to pure tunnelling theory (dashed line) is in qualita- 
tive agreement with experiment. A closer fit is achieved by 
taking into account a low frequency current noise with r.m.s 
value i&I^Y'^'^ ~ 5nA, which broadens the escape curve (thick 
solid line). The effect of a high frequency noise with the same 
r.m.s fluctuations is also plotted (thin solid line). Inset: the 
working point Wp in the bias parameters space. A negative 
flux pulse brings the system close to the critical line /c($(,). 



IV. EXPERIMENTAL RESULTS 

In this section, we present an experimental charac- 
terization of the detection performances. All the re- 
sults have been obtained at the working point Wp = 
(2.222 ilk, -0.117 $o)- The bias flux is close to the one 
which maximizes the critical current (see inset in Fig. 
|4]) and Wp thus presents a low sensitivity to flux fluctu- 
ations. Using the electrical parameters of the SQUID, 
the characteristics of the potential well are calculated: 
Vp = 8.428 GHz and AU = h x 71.1 GHz. The weU con- 
tains 8 energy levels and the anharmonicity is vqi — vi2 = 
160 MHz. The measuring pulse has a duration At = 4 ns. 
Our setup allows for faster measurement, but this timing 
preserves a trapezoidal pulse shape with a well defined 
plateau. This simple shape makes the application of tun- 
nelling theory straightforward. 



A. Detection of the ground state 

Fig. Slshows a plot of the experimental detection prob- 
ability Pg of the ground state versus the measuring pulse 
amplitude. The measured curve is slightly wider than 
the MQT prediction (Eq. [4|). To quantify this, we in- 
troduce the width A<i> of the transition, defined as the 
flux range between po ^ o.05 and P^" = 0.95. MQT 
predicts A$ = 7.46 x 10"'^ $0 whereas the experimental 
value corresponds to 9.52 x 10^^ $o- As discussed below, 
this broadening is an effect of residual noise in the bias 
parameters. 

We restrict the discussion to adiabatic noise, whose 
frequencies are below the typical transition frequency of 
the SQUID, approximately given by Vp. Such a perturba- 



tion does not induce transitions between levels. Without 
any loss of generality, we can focus on a current noise 
61 . The comparison between the fluctuation frequency 
and the timescale of the experiment plays a crucial role 
on the impact of noise. Two limiting situations can be 
distinguished. Fluctuations much slower than At~^ are 
named low frequency noise. In this case, Fq is constant 
during one measurement pulse, but randomly fluctuates 
from one measurement to the other. Assuming that 5000 
repetitions is enough to explore all realizations of 51 (er- 
godicity hypothesis), the actual measured probability is 
the average 

{Pe)if = l-(exp[-F„(4-H5/,$6-l-$™)At]), (6) 

where the brackets denote an averaging over the differ- 
ent values of 61 . As illustrated in Fig. |4j this averaging 
broadens escapes curves. On the other hand, high fre- 
quency noise, with fluctuations much faster than At~^, 
induce random fluctuations of escape rate during one 
measuring pulse. Following Ref. |22l |. the correspond- 
ing high frequency averaging reads 

{P:)hi - 1 - cxp [ - (r„(/b + 61, $fc + $™))Ai] . (7) 

At lowest order, it just translates the escape curve (see 
Fig. [4]) and does not modify its shape. 

Thanks to the massive filtering and shielding of the 
experimental setup, noise on bias parameters is mostly 
generated by the close environment of the SQUID which 
is thermalized at 30 mK. The associated thermal fre- 
quency vt = kBT/h = 600 MHz is much lower than the 
typical transition frequency Vp k, ^ GHz. Therefore, the 
noise can be treated as an adiabatic perturbation. 

In Ref. [ll[, escape measurements of the ground state 
induced by slow current pulses exhibited a low frequency 
flux noise for the same experimental setup and SQUID. 
Its r.m.s. value ^i5<I'^)^''^ = 5.5 x 10~*$0j corresponding 
to noise frequencies ranging from 100 mHz to 20 kHz, is 
largely insufficient to account for the broadening of the 
width. In the following, this flux noise is neglected. 

The close electrical environment of the SQUID also 
generates Nyquist current fluctuations whose spectrum 
density is easily derived from the quantum fluctuation- 
dissipation theorem. Integration between 100 mHz, the 
inverse of a 5000 repetition sequence duration, and 
Ai"i = 250 MHz gives {6PY''^ = 5.5 nA. The best 
flt to experimental data, presented in Fig. HI is ob- 
tained when a low frequency Gaussian current noise with 
((5/^)"'^/^ = 5 nA is taken into account. This value is in 
good agreement with the previous calculation and the 
whole analysis demonstrates we control the measurement 
of the ground state on the nanosecond scale. 

B. Detection in the two level limit 

In the two level limit, the escape probability takes the 
simple form: 

P. = (Pr)+[(P,i)-(P°)] xpi. (8) 



6 



The two detection probabilities (Pe) and (Pg) depend 
on They are calculated taking into account the low 
frequency current noise previously discussed (for clarity, 
the index is implicit). This noise slightly broadens 
escape curves and thus affects the contrast of detection. 
The expected contrast (P^) — (Pe) is 0.95 (compared to 
the 0.98 predicted by the pure MQT regime). We note 
that high frequency noise does not change the separation 
between escape curves of states |0) and |1) and leaves the 
contrast of detection unchanged. 

In the following, |1) is first populated with a low power 
microwave pulse whose 300 ns duration is sufficient to 
reach a stationary state where pi depends only on h' and 
a^w The inset of Fig. [5] shows the resonance curve 
associated with the |0) |1) transition for a microwave 
amplitude a^^, = 16.4MHz. A fit to a Gaussian lineshape 
leads to i^oi = 8.283 GHz and a full width at half max- 
imum equal to 110 MHz. The location of the |1) |2) 
transition, with resonant frequency = vqi — 160 MHz, 
is indicated by the vertical dashed line to the left of foi 
in the inset. This transition is not excited for such a 
microwave amplitude. 

The main plot of Fig. [5] shows the escape probability 
after a microwave pulse tuned to vqi , versus the measure- 
ment amplitude The experimental detection proba- 
bility of the ground state, already discussed, is also plot- 
ted. The contrast of the detection is directly proportional 
to the difference of these two curves. To help interpre- 
tation, the calculated detection probability of |1) is also 
shown. For small |Q) and |1) remain trapped within 
the well. When is increased, a clear bulge appears: it 
corresponds to the selective escape of state 1 1) . At higher 
$m, the potential barrier is so small that both levels es- 
cape. A fit to Eq. [8] quantitatively reproduces the data 
and gives pi = 0.18. The experimental optimization of 
the contrast leads to the amplitude — —0.0754 $o 
which corresponds to a ground state detection probabil- 
ity of about 0.03, a value consistent with theory. In the 
next paragraph, the measurement amplitude is tuned on 
this optimal value. 

To quantify the detection contrast, we increased the 
/iw amplitude to fully saturate the |0) — ^ |1) transition. 
The ampitude of the resonance peak versus the /iw 
amplitude is plotted in Fig. [S] The experimental 
data display two regimes separated by a saturation 
/.tw amplitude about 15 MHz. The low /.tw amplitude 
regime is characterized by a square law dependence 
(solid line). Well above 15 MHz, the resonance ampli- 
tude is nearly independent on the /iw amplitude : the 
|0) |1) transition is saturated. We now focuss on 
the measurement corresponding to a^^, — 50.5 MHz, 
with an associated resonance amplitude equal to 0.30. 
Unfortunately, for such a high power, the /iw excitation 
is less selective and a slight population of |2) occurs. 
A rough estimation using multilevel Rabi theory (T^ 
gives for this contamination p2 — 0.06. The corrected 
experimental population of |1) is thus 0.30 — 0.06 = 0.24. 




FIG. 5: Detection of level 1 1) partially populated with a a^^ = 
16.4MHz amplitude microwave pulse (□). The solid line is a fit 
with pi as the only free parameter, giving pi = 0.18. The ex- 
perimental detection probability of the ground state (o), and 
the calculated detection probability of |1) (dashed line) are 
also plotted. The optimal detection contrast, pointed to by 
an arrow, is obtained for ~ —0.0754 $0; this setting cor- 
responds to Pe ~ 0.03. Inset: spectroscopy of the |0) |1) 
transition for the same /.iw amplitude. The circles are exper- 
imental data and the continuous line is a fit to a Gaussian 
lineshape. The resonance peak is centred on z^oi = 8.283 GHz 
with a 110 MHz full width at half maximum. 




\m amplitude, a^^ (MHz) 

FIG. 6: Amplitude of the resonance peak, such as shown in 
the inset of Fig. [5] versus the /^w amplitude plotted with a 
double logarithmic scale. The measurement amplitude has 
been tuned to maximize the detection contrast. Well below 
a saturation fj,w amplitude about 15 MHz, the experimental 
data (■) display a square law dependence (solid line). In 
the high /xw amplitude regime the |0) |1) transition is 
saturated. 



The theoretical value is 0.50 - 0.06 = 0.44. The ratio 
of these two quantities, equal to 0.54, is the optimal 
contrast of the detection. It is significantly smaller than 
the expected value of 0.95. To understand the reason for 
this contrast loss, we investigate in the next paragraph 
the effects of relaxation during the measurement flux 
pulse. 

Relaxation can occur at three different stages of the 
experiment: between the /.tw pulse and the flux ramping 
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of the measurement pulse (1), during the flux ramping 
(2) and during the measurement stage (3). 

1. At the working point Wp, the energy relaxation of 
state |1) has been experimentally investigated using 
Ti NMR-like measurements. |1) is first selectively 
populated with a 300 ns microwave pulse tuned at 
Vol, with an amplitude a^^, = 16.4 MHz. Then the 
state occupancy is measured with increasing time 
delay after the end of the fxw pulse. The escape 
probability follows an exponential relaxation with 
a rate l/r)^ = 60 ns (not presented here, see Ref. 
[m for a similar curve). Here the delay between 
flux ramping and /iw is kept as short as 1 ns and 
the resulting depolarization is 1.6 %. 

2. During the adiabatic flux ramp, vqi decreases from 
8.28 GHz to 5.42 GHz. During this 2.8 GHz fre- 
quency sweep, the SQUID can couple to spurious 
resonators located in the dielectric tunnel barrier 
of the junctions, as demonstrated in Ref. 14, 2^. 

3. Finally, at the top of the flux pulse, relaxation com- 
petes with tunnelling, yielding a contrast loss given 
by the ratio Tj^/Ti of these two processe rates. A 
fast measurement implies high tunnelling rate and 
thus low depolarization. Indeed, at the measuring 
point, Fi is typically 2.5 GHz. An extrapolation 
of the results presented in Ref. [llj to the transi- 
tion frequency vqi = 5.42 GHz suggests a relaxation 
rate of about 33 MHz. Therefore, the contrast loss 
is very small, about 1.3 %. 

In the following, we neglect depolarization before and 
after the flux ramping. We assume that depolarization 
during the flux ramping is entirely responsible for the ob- 
served contrast loss. We describe this process with the 
help of an effective relaxation rate Feff from |0) to |1) 
which acts during the flux ramping time r = 1.6 ns. The 
population pi (r) just after the ramp is related to the pop- 
ulation pi after the /iw pulse by Pi{t) — p\ exp(— Fofrr). 
Feff, the only free parameter of this phenomenological 
model, is deduced from the experimental detection con- 
trast: 0.95 X exp(-Feffr) = 0.54 gives Foff = 350 MHz. 
In the next section, we show that this model can be ex- 
tended to account for measurement performance when 
more than 2 levels are populated. We will distinguish 
between the set of occupancies just after /iw and 

{p„(t)}, the level populations just after the end of the 
flux ramping. 



C. Multilevel detection 

To investigate the measurement performances in the 
multilevel regime, the states |1), |2) and |3) are populated 
with a Rabi pulse. Fig. [7]^ a) presents a Rabi oscillation 
induced by microwaves tuned on i/qi with an amplitude 
ci/iiu = 164.4 MHz similar to the anharmonicity of the 



(a) 0.8 




state state 



FIG. 7: (a) Experimental Rabi oscillation for a^,„ = 
164.4 MHz. (b) Detection of excited states populated by a 
Rabi pulse with duration = 3.8 ns. The solid line is a 
theoretical adjustment of experimental data (□). For compar- 
ison, the experimental detection probability of |0) is plotted 
(o). The calculated detection efficiencies of pure states |1) to 
j3) are also plotted as dashed lines, (c) Calculated population 
distribution after the end of /^w excitation, neglecting deco- 
herence. During the flux ramping, this distribution relaxes to 
the calculated configuration plotted as solid bars in (d). The 
diamonds are the experimental populations extracted from 
the fit of curve (b). 



potential well. The oscillation is damped in time with 
the characteristic time T2 = 19 ns. 

A microwave duration t^w = 3.8 ns maximizes Pe- In 
the two level limit, this setting corresponds to a Rabi 
TT pulse, which coherently brings the system from |0) to 
|1). With the help of the model developed in Ref. |1^ . 
the population distribution {p^} in the present case is 
calculated (see Fig. [ZKc)). Particularly, the population 
of 1 2) nearly reaches 0.60 and the higher energy state |3) 
is occupied at about 0.10. Fig. [7]^b) shows the evolu- 
tion of the escape probability versus To clarify the 
interpretation, the calculated escape curves of levels |0) 
to 1 3) are also plotted. Surprisingly, there is practically 
no detectable escape of |3). The flrst noticeable bulge 
corresponds to the selective escape of |2) and the second 
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FIG. 8: (a) Measured escape probability after a Rabi tt pulse 
versus microwave amplitude. The measurement pulse ampli- 
tude has been set to maximize the detection contrast between 
|0) and The solid line with solid circles is a calculation 
with no free parameter (details in the text) consistent with 
experimental data (♦). (b) Set of calculated populations after 
the microwave pulse (without decoherence). 



one to the escape of |2) and |1). For higher measure- 
ment amplitudes, aU the states escape. The experimen- 
tal occupancies {pn{T)} just after the flux ramping can 
be extracted from a fit to the following formula: 



Pe = (P^) + E [ (Pe) (Pe) ] >< P« W ' (9) 
■11 = 1 



where the escape probabilities (Pg") depends on 
From the fit, plotted as a solid curve in Fig. [7{b), we 
obtain pi{t) = 0.44, P2{t) = 0.12 and P3(t) « 0. This 
set of occupancies, plotted as diamonds in Fig. U^d), is 
completely different from the predicted values after the 
fiw Rabi pulse. Such a striking discrepancy can be ex- 
plained by a relaxation cascade during the flux ramping. 
In the following, we expand the relaxation model pro- 
posed in the previous section to a 4 level system. 

We keep for |1) the relaxation rate Foff . As n increases, 
the relaxation rate of level \n) dramatically increases. 
The relaxation rate associated with the random transi- 
tion |n) to |n — 1) is proportional to |(n|(^|n — l)| ~^ = n. 
If we neglect other relaxation channels, such as \n) — > 
\n — 2), with a relaxation rate proportional to the anhar- 
monicity coefficient a, the relaxation rate of the nth level 
is nTcff. Thus, populations evolve in time according to 
the closed set of master equations: 





/P3\ 




( -SFeff 








d 


P2 




+3Fcff — 2roff 








dt 


Pi 




+2reff 


— rcff 







\pj 




V 


+roff 










I- 










\Pq / 



which is easily solved analytically. The populations 
{Pn{T)} are obtained starting from the initial popula- 
tions after /iw excitation and letting relaxation act 
during the ramping time r. The result, plotted as solid 
bars in Fig. [Tjd) is in relatively good agreement with 
experimental values. However, our model still overesti- 
mates P2{t)- This may be due to the neglected |2) |0) 
relaxation channel or to decoherence. Nevertheless, the 
essential of the physics is captured by this simple model. 

Fig. IHa) presents the measured escape probability af- 
ter a Rabi tt pulse with increasing microwave amplitudes. 
In this experiment, the measurement pulse amplitude has 
been set to maximize the detection contrast between |0) 
and |1). Our phenomenological model can be used to cal- 
culate the global detection efficiencies of each level and 
to predict the final escape probability given an arbitrary 
initial population distribution {pj^} before the measure- 
ment pulse: 



P, = 0.03 



0.54p? 



0.79p^ 



0.90p°. 



(11) 



(10) 



Without any free parameter, this expression describes 
very well the measured escape probability obtained with 
high /iw amplitude. Relaxation explains why Pg does not 
saturate to exactly 1 when 3 or 4 levels are populated 
after the end of the Rabi tt pulse. We also remark that, 
though levels |2) and |3) undergo fast relaxation during 
the flux ramping, their detection probability is not so 
different from the pure tunnelling prediction given in the 
first part of the article. In fact, to give a crude picture, 
the populations relax towards |1) where they accumulate 
and can be detected. 

The discrepancy which appears for small a^u, is due 
to decoherence: low amplitude ^w implies a slow Rabi 
pulse and bigger errors due to incoherent processes. In 
a two level system, decoherence during a tt pulse leads 
to a population reduction ^ [l — exp(— <^iu/T2)] for the 
first excited state. We use this result to roughly esti- 
mate the complex effect of decoherence in the multilevel 
regime. The previous formula leads to the error bars on 
the predicted escape probability, plotted in Fig. [S^a). 



V. CONCLUSION 

In conclusion, we have demonstrated a very fast 
method for detecting the state of a current biased dc 
SQUID. This technique exploits an adiabatic but ultra- 
fast reduction of the barrier height seen by the trapped 
levels, which allows a selective escape of the excited 
states. Practically, the measurement is performed with 
a roughly trapezoidal flux pulse whose plateau dura- 
tion is a few nanoseconds. The detection probability 
of the ground state is well understood using the stan- 
dard macroscopic quantum tunnelling formula. We re- 
port a slight broadening of the escape curve due to cur- 
rent noise generated by the close electrical environment 
of the SQUID. This broadening is by far insufficient to 
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account for the observed experimental contrast of detec- 
tion in the two level limit. A detailed study of depolar- 
ization processes reveals that spurious relaxation during 
the flux ramping dominates depolarization. The simple 
model that we have developed explains the observed con- 
trast of multilevel Rabi oscillations, when the measuring 
pulse amplitude is set to maximise the detection contrast 
between |0) and Measurement then gives the pop- 
ulation of excited states with a good fidelity. Detection 
contrast between \n) and |n -I- 1) with n > 1, however, 
is strongly lowered by relaxation. Depolarization during 
the flux ramping thus appears as a clear limitation to the 
experimental determination of the levels occupancies in 
the multilevel regime. In phase Josephson circuits, this 
relaxation is attributed to spurious two level fluctuators, 
localized inside the dielectric of the Josephson junctions. 
The remarkable results obtained in Ref. [l^l allows for 
optimism. A new design of the dc SQUID, combined 
with a new dielectric for the junction should lead to a 
significant improvement of the detection contrast. 

We thank W. Guichard, F. W. J. Hekking, L. P. Levy, 
and A. Zazunov for fruitfull discussions. This work was 
supported by two ACI programs, by the EuroSQUP 
project and by the Institut de Physique de la Matiere 
Condensee. 



Appendix 1: from 2D to ID potential 

This analytical calculation is strongly inspired by Ref. 
[l7| , in the general case where the inductance asymmetry 
coefficient rj is finite. Let us consider the two dimensional 
potential: 

U{x,y) = 2E j[- sx - cos X cos y-s'qy + h{y~yBf], (12) 

where s = Ib/2Io and ys = 7r($b/$o) are the reduced 
bias parameters. For the bias points used in the arti- 
cle, the potential surface displays equivalent local min- 
ima separated by saddle points. The critical current of 
the SQUID is geometrically defined as the bias current Ic 
where the minimum vanishes and merges with adjacent 
saddle point at the coordinate {xc,yc)- Experiments are 
performed close to the critical current, therefore a Taylor 
expansion of the potential in the vicinity of the critical 
point {sc,Xc,yc) is justified. Following [l7|, we introduce 
Ax = X — Xc et Ay = y ^ y^. Limiting calculation to the 
third order one obtains: 

U/i2Ej) = (s,-s)(Ax + 77Ay) 

+ ^ cos Xc COS yc [Ax^ + Ay^) 

— sin Xc sin ycAxAy + bAy^ (13) 

— i sin Xc cos yc (Ax^ -I- 3Ay^) Ax 

— ^ sin yc cos Xc ( Ay^ -I- 3Aa;^) Ay. 

We are interested in the dynamics along the escape di- 
rection which corresponds to the direction of minimum 
curvature at the critical point. It makes an angle 9 with 



the X axis: 



tan0 — — 



dlU{xc,yc) _ cos Xc cos yc 
dxdyU{xc,yc) sin Xc sin yc 



(14) 



Let (f) be the phase along the escape direction. Looking 
for extrema of the potential along this direction, one finds 
a minimum located at (/> = and a saddle point. Along 
the escape direction, the anharmonic potential reads: 



6AUY 



(15) 



It is fully determined by the plasma frequency Vp 
u!p/{2tt) and the barrier height AU: 



AU ^ (—Y^^ (cose + Tjsmer'' AW (lb, Ic) (16) 

Vug/ 



'^p — — - j (cos tf + rj sm ( 



)3/2 ^^JJ 

^l/4 



iy'p'ih,Ic,2Co). (17) 



The prefactor — sin Xc cos yc cos 0(1 + 2sin^6') + 
sinycC0SXcSin6'(l -I- 2cos^0) is the third derivative of 
the potential along the escape direction. AU^'' and 
are the well known barrier height and plasma frequency 
for a single Josephson junction with critical current /c, 
capacitance 2Co, biased with a current 7^: 



AU"{hJc)^{^^h){l~hlh 

V2 1 Ic \l/2 



3/2 



(18) 



Appendix 2: coupling to a flux signal 

In the general small flux perturbation 5^{t) 

added to the bias flux introduces a perturbation term: 



"I>n 1 

W{t) = - — 5^{t) X y. 



(20) 



The component of this operator on the escape direction 
reads: 



sin 6* 1 / 2/i \i/2_ ^ ^ 
.{t)^^ — —[—-) Sm^V- (21) 



Assuming a monochromatic perturbation 5^{t) = 
^ fj_w cos(2'Kiyt) , one obtains the perturbation term: 



Wip{t) = —hafj_u,V2 cos{2Tn't)ip, 



(22) 



where a^^, is the 2-level Rabi frequency for a tuned exci- 
tation {v = i>Qi): 



sin 61 1 1 



2n Ls yJC^hvp 



(23) 
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